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■ Abstract 

^ ' As a contribution to the inverse scattering problem for classical chaotic systems, we 

^ ■ show that we can select sequences of intervals of continuity, each of which yields the 

i information about period, eigenvalue and symmetry of one unstable periodic orbit. 
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I In the framework of the study of the inverse scattering problem for classical 

Q I chaotic scattering in two dimensions, the attention has been focussed on the 

topology of the chaotic saddle [1,2] and on the definition of an appropriate 
O ■ partition [3]; thermodynamic quantities were also discussed marginally [1]. In 

>■ ! none of these approaches the specific properties of the unstable periodic orbits 

^ I were explicitly searched for. Yet knowledge of these, goes a long way towards 

^ I understanding the chaotic saddle. Usually the shortest orbits overshadow the 

chaotic set or, in other words, they form the skeleton of the globally unstable 
component of the invariant set. As far as semi-classical approximations are 
concerned, these orbits form the backbone of all considerations that lead to 
trace formulae [4-6]. The most useful informations are their periods and Lya- 
punov exponents. These can be used to determine the hierarchical order in 
the pattern of intervals of continuity (henceforth abbreviated to IOC) in the 
scattering functions, but at the same time it is this very pattern, which allows 
us to learn something about the periodic orbits from scattering functions. 

In this article we shall focus on the latter part, i.e. on obtaining the periods 
and the Lyapunov exponents of unstable periodic orbits from regular patterns 
of the IOC of scattering functions. The latter we define as a function of the 
points on a line in the space of initial conditions, that gives some property of 
the scattering process such as e.g. time delay or scattering angle (cf fig. 1). 
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After recalling some known results [7] concerning the external periodic orbits 
we present two methods to achieve our goal for periodic orbits that are both 
fairly short and not too unstable. We illustrate these methods by scattering 
a charged particle off a magnetic dipole. The reader may want to refer to the 
figures of the example while reading the description of the general method. 

The trajectories belonging to one IOC of a scattering function all leave the 
interaction region crossing the same external orbit i.e. leaving through the 
same saddle or "doorway" . When the initial condition approaches one of the 
boundary points of the IOC from the inside, the scattering trajectory executes 
an increasing number of revolutions along the saddle orbit before it finally 
leaves the interaction region. For each revolution we expect an oscillation of 
the scattering angle (cf fig. 2). The frequency of this oscillations increases with 
the diverging time delay. An asymptotic observer can determine a sequence 
of initial conditions from within one and the same IOC, which are separated 
by a full oscillation of the scattering angle. This sequence converges to the 
boundary of the IOC in such a way, that, in the limit, the distance from 
the boundary point decreases by a factor determined by the eigenvalue of 
the saddle orbit. In the same limit the time delay increases by the period of 
this orbit, when going from one member of this sequence to the next. Prom 
these two quantities the Lyapunov exponent can be calculated. If there are 
several external orbits we can investigate each of them separately, choosing 
an appropriate IOC, where the outgoing scattering trajectories cross exactly 
this orbit before finally leaving the interaction region. Information about any 
of the internal orbits cannot be obtained in this way. 

Based on this result we propose two methods to obtain both periods and 
Lyapunov exponents for inner periodic orbits. The first will use the time de- 
lay function only whereas the second exploits in addition the other scattering 
functions. To achieve this we analyze how the inner orbits influence the hier- 
archical pattern of IOC. 

The basic idea for our method is provided by the following property of chaotic 
scattering: To each periodic orbit a, which can be reached from the outside, 
there exist infinitely many sequences of IOC (cf fig. 1) of the scattering func- 
tions. The common feature of these sequences is the following: As we move 
along the sequence, sooner or later one step implies precisely one revolution 
of the trajectory close to the periodic orbit a. The central problem will be to 
identify a sequence which is clearly associated with one periodic orbit. If we 
knew such sequences in advance, then it would be easy to find the period and 
the Lyapunov exponent of the periodic orbit a. 

We label the IOC forming the selected sequence by {/^'^'^jfegN where k labels 
the order of the elements of the sequence, L[l'^'>) indicates the length of the 
IOC and t{l'^'>) a typical time delay for this IOC which we obtain by choos- 
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ing a representative trajectory for each IOC. It is convenient to choose the 
representatives to have the absolute minimal time delay within an IOC. 

At large time delays two consecutive representatives differ only in one addi- 
tional revolution along the unstable periodic orbit. In order to extract peri- 
ods and Lyapunov exponents of unstable periodic orbits from the time delay 
function, wc need to know the representative time delays and lengths of as 
many IOC as possible. We determine the ratio A^''^ = L(/^''))/L(/^''+^)) of the 
lengths of consecutive pairs of IOC and the difference of their representative 
time delays fj*^) = ^(4''^^^) - ^(4*^^)- We also define 

If the sequence was chosen appropriately, T^''^ will converge to the period T^. 
of the orbit a, X^^^ to its eigenvalue Ao- and A^.*^) to its Lyapunov exponent A^-. 

As mentioned above the choice of the sequence is critical and difficult. We 
therefore proceed to show how this choice can be partially avoided simulta- 
neously with the approximative evaluation of the quantities above. If we plot 
either A^*^^ or A^*^) against T^*^^ we should find an accumulation point near the 
correct value of {Xa,Ta-) or (Aa^T^) respectively. Note that each of the infinitely 
many sequences associated to one periodic orbit will approach the same ac- 
cumulation point. Furthermore points (Ao-',To-') or (Ao-',To-') erroneously taken 
because they belong to a sequence for a different periodic orbit will have dif- 
ferent accumulation points, except if both period and eigenvalue for the two 
orbits coincide. If the error is such that the pairs are obtained from IOC of 
sequences that belong to different periodic orbits no accumulation points are 
expected. It is thus tempting to perform this plot using points generated for 
all pairs of IOC available. 

In principle we should then find accumulation points for many periodic orbits. 
This will not be implemented easily because if we have sufficient IOC to ex- 
pect reasonable convergence along one sequence, the plot would have so many 
points that it would be very difficult to identify these accumulation points. 
Even worse, the discrete numerics may simulate structures that do not exist. In 
order to identify the accumulation points, only pairs of IOC that might belong 
to the same sequence have to be taken into account. Appropriate strategies to 
achieve this goal must be developed. An inspection of the hierarchical struc- 
ture corresponding to the time delay function shows that successive IOC have 
to be connected by pieces of equal or higher hierarchical order. Therefore only 
pairs of IOC are tested with no intervening IOC possessing lower time delays. 
Further filtering may well be necessary. 

• If we proceed using the time delay function as the only source of information. 
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we can start to look for multiplets of IOC which apparently belong to one 
sequence. Each of its members should have a higher representative time 
delay than its predecessor without intervening IOC with lower time delays. 
If such a multiplet belongs to a single sequence its successive IOC should 
show roughly the same differences in time delay and ratios of successive 
lengths (cf fig. 3). By narrowing the limits for the allowed differences in 
time delay and ratios of the multiplets, improved convergence towards a 
periodic orbit can be promoted. 
• If other scattering functions are also available, pairs of IOC that belong to 
a sequence can be identified by comparing their representative final con- 
ditions. Again only pairs of IOC without intervening IOC with lower rep- 
resentative time delays are taken into account. Since the final conditions 
of representative trajectories converge in the limit of high time delays to a 
point, two successive IOC can be tentatively identified by looking for those 
with close representative final conditions. 

Both methods yield not only the basic periods but also, in a diminished ex- 
tent, integer multiples thereof. 

To illustrate these methods we use the scattering of a charged point parti- 
cle off a magnetic dipole. This model is also known as the Stormer problem 
[8]. After eliminating the rotational degree of freedom and a suitable rescaling 
in cylindrical coordinates one is left with a Hamiltonian of the form: 



Here {pp,Pz) are the momenta conjugate to {p,z). It is well known that the 
resulting phase fiow represents a chaotic scattering system [9,10]. An analyti- 
cal proof for its non-integrability is given in [1 1] . Every scattering trajectory 
is uniquely determined in the asymptotic region by the impact parameter b^^ 
and the angle a^^ between its velocity vector and the negative p-axis. A de- 
tailed description of the chaotic dynamics of this system for the energy range 
[0.031, 0.081] can be found in [10]. The two fundamental periodic orbits of the 
binary horseshoe will be called A and C. The outer orbit C is accessible from 
the outside for all energies whereas the inner one A may be at the center of 
an elliptic island for some energies. 

For illustration we will consider the hyperbolic case at an energy E = 0.062 
and show how to extract information on A and other inner periodic orbits. 
Note that hyperbolicity is no prerequisite for the methods presented here. For 
the integration of the equations of motion a Bulirsch-Stoer algorithm [12] has 
been used. The simulation takes place in a disk of radius 1000 around the 
origin of the (p, 2;)-plane. 

A part of a typical sequence whose trajectories come close to A is shown as 
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IOC in a time delay function in fig. 1. In fig. 2 the final angles ccout oi the 
trajectories from the IOC {i^^ '*}i<jk<5 are plotted against their corresponding 
time delay values. 

In this plot several essential features can be seen. 

(1) As the time delay increases the period of the curves converge towards the 
period of the outer periodic orbit C. This is a direct consequence of the 
additional revolutions near C with initial conditions near the boundary 
of an IOC. 

(2) The differences of the representative (minimal) time delays of con- 
secutive IOC converges towards half the period of the inner periodic orbit 
A. 

(3) The lower end of the curves of every second IOC point in the same di- 
rection. This effect has its origin in the common reflection symmetry on 
the p-axis of A and the Hamiltonian. If A and the Hamiltonian had no 
common symmetry 

t(/?^)=^(t(/r^)+t(/r^)) (3) 

would not be fulfilled. In such there would be two sequences whose 

members were separated by the full period of the corresponding periodic 
orbit. This shows that by using only the time delay function as a source of 
information it might be unclear whether we measure the whole period of 
a periodic orbit without symmetries or its n-th part of an periodic orbit 
with a cyclic symmetry group of order n. Using additional information 
provided by other scattering functions as ctout and bout can give helpful 
hints to the answer of this question. 

In order to resolve about 608 IOC we scatter 10^ trajectories with = 0° 
and 6in evenly distributed on the interval [1.3049,1.3065]. If we apply the 
filter that respects the hierarchical structure and look for triples of IOC 
(/^^~^), 4''+^)) with the property 



f{k-i) __ f(k) 

a cr 



<0.02 
<0.02 



(4) 



we end up with 64 of such triples. Each triple of IOC yields two points 
(AW^fW) in fig. 3. Not all 128 points are located in the area shown in this 
figure but the chosen resolution of 10^ points was too poor for a meaningful 
prediction of orbits with symmetry reduced periods longer than T = 15. The 
more dots are clustered in a group the likelier it is to find a periodic orbit with 
roughly these values for its symmetry reduced period and Lyapunov exponent. 
In fig. 3, apart from the times for one revolution around a symmetry reduced 
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periodic orbit (A^/^, C^/^) integer multiples of these times are also recogniz- 
able {A, C). Among the periodic orbits of our system shown in fig. 4 only three 
orbits {A,C and /) remain invariant under a reflection on the p-axis. The rest 
of the orbits either change their direction of rotation {B,D,F,G and H) or 
are mapped onto their mirror image {E). Thus for orbits A,C and / the half 
periods show up besides the full ones in fig. 3. 

Alternatively we may use pairs with similar final conditions retaining the re- 
striction derived from the hierarchical structure of the time delay function. 
We accept pairs of IOC whose representative trajectories differ in outgoing 
angle by less than 0.1" and in the outgoing impact parameter by less than 
0.01 in our scale. We plot the 115 obtained points in fig. 5 similar to fig. 3. 
Now half periods will not appear because we did not symmetrize the outgo- 
ing asymptotic parameter. Comparing the two figures we see that the second 
method allows to identify more periods. While in fig. 3 accumulation points 
are not very obvious in fig. 5 they appear rather clearly in a few instances. 
Where they appear the Lyapunov exponent indeed seems to converge toward 
the exact value. Note the difference in scales for time delay and exponents. 
Considering that the Lyapunov exponent is an average, once identified, the 
points belonging to one period can be averaged. Table 1 shows that for the 
points in fig. 5 the average matches the exact value quite well. 

Summarising we have presented a new approach to the inverse scattering prob- 
lem for chaotic Hamiltonian systems. In distinction to earlier work [1-3] we 
do not require the system to have an internal or external clock. Furthermore 
we are able to extract more detailed information about the properties of the 
most important periodic orbits. The possibilities inherent to this approach are 
much larger than what was presented in this letter. More extensive research 
on the filtering techniques must be done and the optimal approach may well 
depend on the problem. Also much larger numbers of intervals of continuity 
must be generated to be able to analyze the structure of the (A, T')-plane. Fi- 
nally the possibility to obtain quite easily the periods of the shortest periodic 
orbits may be used to provide an inner clock that always exist for the method 
proposed in [1]. 
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Fig. 3. Periods and Lyapunov exponents of periodic orbits evaluated by only using the time delay function as described in the text. Exact 
values are marked by open circles and capital letters. The letter's upper indices display what fraction of the full period were measured. 
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Fig. 5. Periods and Lyapunov exponents of periodic orbits evaluated by using the time delay, CKout and 6out scattering functions. Exact 
values are marked by open circles and capital letters. The letter's upper indices display what fraction of the full period were measured. 
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Table 1 

Measured Periods and eigenvalues of periodic orbits. Sample region was E = 0.062, a-m = 0, bin € [1.3049, 1.3065]. RMS errors. 



